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DIFFUSION  OF  MAGNETIC  FLUX  INTO  A CYLINDRICAL  LINER  IN  THE 
LIMIT  OF  HIGH  MAGNETIC  REYNOLDS  NUMBER 


I.  Introduction 

Over  the  past  two  decades,  a number  of  laboratories  have  carried 
out  magnetic  flux- compress ion  experiments  using  imploding  conducting 
liners.1-7  The  liners  are  cylindrical  and  surround  flux  lines  parallel 
to  the  axis  of  symmetry.  The  axial  flux  may  be  accompanied  by  a gas  or 
plasma  fill;  in  some  experiments  the  latter  constitutes  most  of  the 
payload  (i.e.,  supplies  most  of  the  pressure  opposing  the  liner 
implosion),  while  the  magnetic  field  serves  only  to  reduce  thermal 
losses  or  as  a diagnostic. 

The  liners  have  been  formed  of  metal  in  either  solid  or  liquid 
state.  Chemical  explosives,1*2  electromagnetic  acceleration  in  the 
form  of  9-pinch3,s  and  Z-pinch4“s  drives,  and  hydraulic  compression7 
have  been  used  to  implode  the  liner.  With  these  three  driving 
agencies  are  associated  typical  liner  velocities  v of  order 
106  cm/sec,  105  cm/sec  and  104  cm/sec  ("fast",  "intermediate",  and 
"slow"  liners  ),  respectively.  (These  velocities  scale  roughly  as 
v^  ~ (E^/p^)-1-/2,  where  E^  is  the  energy  density  in  the  drive  and 
P is  the  liner  mass  density.) 

A limit  on  the  value  of  the  field  strength  B which  can  be 

obtained  is  imposed  by  diffusion  of  flux  from  the  compression 

volume,  owing  to  the  finite  conductivity  of  the  liner.  If  the  liner 

is  given  an  energy  per  unit  length  W /I  and  all  of  this  goes  into 

Li 

flux  compression,  then 

Wl/4  » nr2B2/2|A0, 

Note:  Manuscript  submitted  May  26,  1977. 
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where  r is  the  final  compression  radius.  A skin  depth  A is  associated 
with  the  diffusion  of  flux  lines  into  the  liner  over  a time  t 
according  to  A2  = K^t,  where  kq=  l/^M^  is  the  dlffusivity,  is  the 
electrical  conductivity,  is  the  (free  space)  magnetic  perme- 

ability , and  t » r/vL  is  the  time  scale  of  the  liner  motion.  The 
condition  that  the  flux  which  leaks  out  be  less  then  that  inside  the 
compression  volume  is  roughly 


2nrA 

nr2- 


< 1. 


(2) 


This  is  simply  a requirement  that  the  magnetic  Reynolds  number  Rm, 
defined  by 


Re 


(3) 


be  greater  than  unity.  For  a conductivity  equal  to  that  of  liquid 
copper,  Eqs.  (l)  and  (2)  imply 


B < 17  (WL/f)l/avL, 


(*0 


where  B is  in  MG,  W /f  in  Mj/m,  and  v is  mm/M'Sec. 

L L 

In  the  present  paper  we  describe  several  calculations  intended  to 
elucidate  the  roles  of  (i)  convergence  (thickening)  of  a cylindrical 
liner  and  (ii)  the  detailed  time  dependence  of  liner  trajectories  in 
the  diffusion  of  magnetic  flux  into  a metal  liner.  These  calculations 
have  been  carried  out  in  support  of  the  NRL  Linus7'8  experiments,  in 
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which  a rotating  liquid  metal  liner  is  imploded  hydraulically  onto  a 
mixture  of  plasma  and  magnetic  flux.  Ideally,  nearly  all  of  the 
compressional  energy  should  go  into  the  plasma,  thereby  heating 
it  to  kilovolt  temperature  and  initiating  thermonuclear  reactions, 
and  the  magnetic  field  should  be  restricted  to  a thin  buffer  layer 
between  the  plasma  and  the  liner.  Maintenance  of  this  layer  for 
the  duration  of  the  burn  is,  of  course,  essential  to  containment. 
Even  if  only  a small  portion  diffuses  into  the  liner,  the  Joule  heat 
produced  by  the  penetration  of  current  into  the  liner  can  vaporize 
the  latter  and  release  high  Z impurities  into  the  plasma.  There- 
fore, a reliable  quantitative  theory  of  the  diffusion  is  necessary. 

In  current  and  near-term  experiments,  however,  much  or  all  of 
the  "payload"  consists  of  magnetic  flux.  In  this  paper  we  therefore 
simplify  our  problem  by  ignoring  the  plasma  entirely.  We  also  ignore 
field  line  fringing  and  curvature  so  as  to  make  the  problem  one- 
dimensional. (This  is  strictly  correct  at  the  midplane  of  the  experi- 
mental device,  if  nowhere  else.) 

Within  the  limitations  thus  invoked,  our  objective  is  to  cal- 
culate in  detail  how  much  flux  leaks  out  of  the  compression  region, 
retaining  realistic  conductivities,  liner  trajectories,  and  the 
associated  liner  convergence.  Previous  work9  has,  for  the  most  part, 
been  restricted  to  slab  geometry,  using  a vacuum  magnetic  field 
prescribed  on  the  basis  of  analytical  convenience  rather  than  con- 
sistency with  realistic  trajectories. 


I 
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Below  in  Section  II  we  develop  a formalism  applicable  to  either 
slab  or  cylindrical  liners.  A detailed  treatment  of  the  former  is 
presented  in  the  Appendix  to  facilitate  reference  and  comparison 
with  the  cylindrical  case.  In  Section  III  we  specialize  to  the  limit 
of  greatest  interest,  the  weak-diffusion  (high  magnetic  Reynolds 
number)  limit,  and  solve  the  resulting  boundary  layer  problem.  The 
calculation  of  total  diffused  flux,  magnetic  energy  in  the  liner  and 
Joule  heat  is  reduced  to  a series  of  quadratures  with  respect  to  the 
trajectory.  The  diffused  flux  per  unit  surface  area  is  precisely 
the  same  as  in  the  slab  problem.  In  other  words,  the  flux  diffusion 
rate  is  proportional  to  the  instantaneous  radius.  Therefore  the 
total  diffused  flux  is  proportional  to  an  average  radius,  whose  value 
depends  on  the  details  of  the  liner  trajectory.  Section  IV  describes 
the  results  of  applying  this  formalism  to  two  prescribed  model  tra- 
jectories, termed  "sinusoidal"  and  "parabolic."  We  conclude  in 
Section  V with  a brief  discussion  of  the  results. 


II.  Derivation  of  the  Formalism 

Consider  a situation  in  which  a cylindrical  vacuum  region  of 
radius  Ri  contains  a uniform  axial  magnetic  field  (see  Fig.  l).  The 
material  surrounding  the  vacuum  (the  liner)  is  assumed  to  be  incompressible, 
extends  to  infinity  and  has  a finite  electrical  conductivity  a.  Rx  may 
be  a function  of  time,  in  which  case  there  is  a radial  flow  field  with- 
in the  medium.  It  is  assumed  that  at  t = 0 the  magnetic  field  vanishes 
outside  the  vacuum  region.  Its  subsequent  evolution  is  determined  by 
Ampere ' s law 

V x B = MJ  , (5) 

Faraday ' s law 

6B 

“St-  = " ~ x £ y (6) 

and  Ohm's  law 

ctCe  + v x b]  = j , (7) 


where  v is  the  velocity  field.  With  the  assumptions  of  one-dimensionality, 
incompressibility  and  straight  field  lines,  Eqs.  (5-7)  yield  the  con- 
vective diffusion  equation 


&B 

TiT  = 


6 

TT 


(nr 


6 -l 


t>B  >. 

-rr>> 


(8) 
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where  6 = 1,2  for  slab,  cylindrical  geometry,  respectively,  and  h is 
the  diffusivity,  defined  by 


K = 1/p-CT  . (9) 

In  general  h depends  on  position  explicitly  or  through  the  temperature 
dependence  of  a,  and  we  can  write 

H = H F . (10) 

o 

Here  k is  a const  " indicative  of  the  scale  of  h and  F is  a 
o 

dimensionless  which  can  depend  on  r and  the  history  of  B. 

Corresp  uo  Eq.  (8)  the  initial  condition  is 

B (r,0  ) = 0,  r > Rx  (0  ) , (ll) 

and  the  boundary  conditions  are 

lim  B ( r , t ) = 0 ; (12) 

-♦  00 

CO 

Bi(t)Ri6(t)  + bj  drr6_1B(r,t)  = Bi(o)Ri(o)  • (15) 

Ri 

Equation  ( 13 ) is  the  expression  of  conservation  of  total  magnetic 
flux.  These  boundary  conditions  correspond  to  an  initially  field- 
free  liner.  The  case  where  there  is  a non-vanishing  field  profile 
within  the  liner  initially  can  also  be  treated,  and  the  results  can 
be  related  to  those  of  the  present  case.  For  large  compression  ratios 
(see  Section  IV ),  the  effect  of  the  initial  profile  is  negligible. 


w 


It  is  convenient  to  scale  B(r,t)  by  the  maximum  value  of  the 
vacuum  field  associated  with  the  implosion,  Bq: 

b = B(r,t)/Bo,  (14) 

and  to  scale  t by  a time  tQ,  characteristic  of  the  liner  motion, 


T = t/tQ  • (15) 

Equation  (8),  which  represents  the  Eulerian  form  of  the  convective- 

diffusion  equation,  is  not  in  the  most  advantageous  form  to  calculate  the 

diffusion  into  a moving  liner.  This  is  because  the  boundary  condition 

of  total  flux  conservation,  Eq.  (13),  is  applied  at  the  moving  boundary, 

Rx(t).  Furthermore  there  are  really  two  time  scales  associated  with 

« 

Eq.  (8):  the  convective  time  scale,  Ri/Ri(t),  and  the  diffusive  time 
scale,  k /Rj.2.  It  is  most  useful  to  cast  Eq.  (8)  in  Lagrangian  form. 

For  the  incompressible  liners  considered  here  the  Lagrangian  variable  is 

5 - (r6  - Ri)/Ro6  . (16) 

where  Rq  is  the  minimum  radius . The  system  of  equations  now  becomes 


6b  - Rm-l  b 

■F’"8"  "51 


F[Rx6  + g]  6 

6b  \ 

6?  | 5 

(17) 

? > o ; 

(18) 

! 0 ; 

(19) 

?-♦<*> 


7 


CO 

Mt)R I6  +J  dSb(?,T)-l, 


(20) 


where  ^(t)  = b(0,T)  = Bi/B^  the  scaled  liner  radius  Ri/RQ  is  denoted 
by  R1#  and  the  magnetic  Reynolds  number  Rm  is  defined  by 


Rm  = R2/52h  t 
O o o 


(21) 


The  slab  case,  shown  in  Fig.  2,  is  treated  in  the  Appendix. 

Equations  (17-20)  with  6=3  (spherical  geometry)  apply  to  such  processes 
as  heat  diffusion  from  a collapsing  reservoir,  but  not  to  magnetic 
flux  diffusion.  Spherical  (or  quasispherical)  flux  surfaces  not  only 
do  not  arise  even  in  three-dimensional  implosions,  they  also  contradict 
the  assumption  of  straight  field  lines  employed  in  the  derivation  of 
Eq.  (8).  Thus  in  the  body  of  this  paper  we  restrict  ourselves  to  the 
cylindrical  case,  6=2,  and  rewrite  the  equations  accordingly.  The 
Lagrangian  spatial  variable  is  then  from  Eq.  (16) 


5 = (r2-Rl2)/Ro2, 


(22) 


and  the  appropriate  form  of  the  equations  is 


6b  _ 6 r 2 . a \ 6b 


6 t 


b(5,0)  = 0, 


*|- [<Sx8  + 5)-g-]; 


§ > 0; 


(23) 

(2k) 


lim  b(5,T  ) = 0 ; 

5 -»  oo 


(25) 


bi(T)Ri2+J  dS  b(S,T)  = 1, 


Rm  = R 2/4  h t . 
o o o 


Some  results  of  an  approximate  nature  can  be  found  when  Rm  = 0 
(static  liner).  Given  the  value  of  b at  r = R1;  the  solution  for 
the  diffusion  of  b outward  from  the  surface  r = Rx  into  a region  where 
b = 0 initially  is  well-known10: 


b(P,T)  = dr'b1(r')Jdu  u e“u  (T'T 


VU)  VpU)  " J0(pu)  Yo(u) 

j2(u  ) + Y 2(U ) 

o o v ' 


where  P = r/Rx . Substitution  into  Eq.  (26)  yields  an  integral  equation 
for  ba. 


bx(T)  + ( | )2f  dr'biCT-T  *)j 


= 1 (29) 


For  small  t bj  1 , and  the  lowest  order  correction  can  be  found  by 

3 

iteration,  using  the  asymptotic  expansion  for  the  Bessel  functions 

J and  Y : 
o o 


bl^)-  1 - ( |)3/TdT^  - 


-U 


* l - 


9 


Analogously,  for  t -*  » the  small-argument  expansion  of  the  Bessel 
functions  leads  in  lowest  order  to 


bi(T)~;r-krF  • (51) 

We  note  that  bz (t  ) -*  0 faster  than  in  the  slab  case  (see  Appendix), 
in  accord  with  the  intuitive  notion  that  radial  efflux  is  less 
constrained  than  is  rectilinear  efflux. 


10 
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HI.  Large  Magnetic  Reynolds  Number 

The  case  of  chief  experimental  interest  is  that  of  large  magnetic 


Reynolds  number, 

Rm  = R 2/4h  t » 1,  (32) 

since  this  is  the  limit  of  efficient  compression.  Now  the  situation 
is  greatly  simplified.  We  anticipate  that  diffusive  effects  should 
be  localized  in  a very  thin  boundary  layer  near  the  liner-vacuum 
interface.  That  is,  the  magnetic  flux  distribution  should  differ 
only  by  a small  perturbation  from  that  for  an  ideal  adiabatically 
compressed  trapped  flux. 

To  establish  this  result  rigorously,  we  proceed  using  the  methods 
of  matched  asymptotic  analysis.  The  relevant  equations  are  Eq.  (23), 
the  Lagrangian  form  of  the  convective -diffusion  equation,  with  initial 
conditionsgiven  by  Eq.  {2k)  and  boundary  conditions  given  by  Eqs.  (25) 
and  (26).  It  is  assumed  in  the  analysis  that  the  form  factor  F remains 
of  order  unity.  We  identify  and  treat  separately  two  regions,  as 
follows. 

Outer  Region 

We  assume  a regular  perturbation  expansion  for  the  scaled  magnetic  field 
strength  b in  powers  of  the  small  parameter  Rm  1^2, 

b(?,T  ;e;Rrn)  = b^°)(?,T;e)  + Rrn-1^2  b^1)  (S,T;e)  (33) 


+ Rm-:!  b(2)  (§,r;e)  + ... 


To  lowest  order  in  Rm"1/2  we  find  from  Eqs.  (23-26) 


b(o)(5,0)  =0,  ? > 0; 
lira  b(°)(0,T)  = 0; 

£ -4  °o 

b(o)(0,T)  = bX(T)  = Rx-2  . 


(34) 

(35) 

(36) 

(37) 


The  solution  of  Eq.  (34)  which  satisfies  conditions  (35)  and  (36)  is 

P 

b(o)(?,r)=0.  (38) 

However,  this  solution  does  not  satisfy  the  second  boundary  condition, 
Eq.  (37).  We  thus  have  the  ingredients  for  a classical  boundary 
layer.  It  is  necessary  to  rescale  the  independent  variable  ? to 
recover  the  missing  boundary  condition  and  resolve  this  layer. 

Inner  Region 

The  rescaling  is  performed  by  defining 

11  = Rm1/2^  . (39) 

When  Eq.  (39)  is  introduced  in  Eqs.  (23-26)  we  obtain 

“It-  = [E(Ri2(t)  + Rm  l/2Tl)  -5^-]  ; (4o) 
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MV>)  = o ; 


lim  b(V)  = 0 ; 


(U 

b(o,T ) + Rm  1^2j'  d”nb(1l,T)  = 1. 


Again  expanding  as  in  Eq.  (33)  yields  in  lowest  order 


yr R1  “5TT  ( 5^1 j ’ 


b(o)n,0)  = 0 ; 


lim  b(°)(T|,T ) = 0; 
T}-*® 


b(o)(0,T)  = bi(T)  = Ri"2  . 

This  is  just  the  condition  for  adiabatic  flux  compression. 
If  we  define  a new  timelike  variable,  the  "age" 


e R^tOdT' 

•'n 


Eq.  (Ml)  simplifies  to  a slab-type  diffusion  equation, 


6b(o)  , 6 6b(°>  , 

"55 5T  <F  TH — >• 


Equation  (49)  embodies  both  the'  cylindrical  geometry  and  the  liner 
motion  in  a relation  which  is  formally  identical  with  that  obtained 


i 
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in  planar  geometry.  Thus  for  Km  » 1,  any  calculation  which  can  be 
done  in  slab  geometry  can  be  carried  over  to  the  cylindrical  case. 

If  at  this  point  we  rule  out  the  possibility  of  thermal,  etc., 
variation  of  o and  set 


F = 1, 


we  can  use  the  results  derived  in  the  Appendix  to  write  down  the  diffused 
magnetic  field. 


B(T],e)  = 


»o  f 

J*  J 


ds  e"s  bx(r  _ T)2As2), 


Tl/a/e 


where  b^.  is  given  by  Eq.  (47).  From  Eq.  (50)  we  find  the  total  diffused 


flux. 


i - 2Tt  I rdrB  ( r,  T ) 


\/tt  B R 2 
o o 


f ' bl 

- de  — 

J ye3 


bi(e') 


the  total  magnetic  energy, 


\-2af 


n3/2  BAa  . rQ  „ b1(e')b1(e<r) 


/-7‘ 

o *o  * 


/ 2 e-e  -e 


!•  ! 


and  the  total  Joule  heat. 
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w 


T 


B 2R  2 
o o 

Rmi72r" 


c 

/ 


de'bi(e') 


i_ 


Bbi/^e"  - 

(e  '-e^2 


(53) 


Equations  (50-53)  correctly  determine  $,  and  WT  for  any  form  of 
the  vacuum  field  Bi(t),  even  if  Eq.  (47),  the  condition  for  adiabatic 
flux  compression,  is  not  satisfied.  In  the  present  instance, 
bi(9)  = [Ri(9^]”2  and  therefore  the  liner  field,  diffused  flux,  and 
magnetic  and  thermal  energy  distributions  are  seen  to  be  given  as 
quadratures  over  the  liner  trajectory,  Ri(T),  in  this  large  Reynolds 
number  limit.  Now  all  that  remains  is  to  supply  the  trajectory  Ri(t). 
In  general  Ri(t)  must  be  obtained  by  considering  all  the  forces  acting 
on  the  liner.  To  emphasize  the  physics  of  the  problem  in  its  most 
simple  form,  we  will  instead  consider  the  radial  trajectory  to  be 
prescribed  and  generate  solutions  consistent  with  this  trajectory. 
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IV. 


Numerical  Results 


/ 


We  wish  to  apply  the  preceding  results  to  trajectories  which  are 
symmetric,  smooth  and  contain  the  compression  ratio  or  (initial  radius 
divided  by  final  radius)  as  a parameter.  Two  examples  are  the 
sinusoidal  trajectory, 

RQ2/R12;t)  = e + (l-e)  sin  (Ttt/2t  ),  (54 ) 


r 


and  the  parabolic  trajectory, 

Ri2  )/Rq2  = 1 + ;e_1  - l)(l-t/tQ)2.  (55) 

In  both  of  these  R’  = R , the  minimum  radius,  when  t = t , and 
_o’  ’o’ 

/ -1/2 

R-i  0-  = e~~'2  R , i.e.,a  = e 
o’  ’ 

But  there  is  an  essential  difference.  For  the  sinusoidal  model, 
RitQ/Ri  =»  - tt/4-s  » 1 at  t =0,  while  in  the  parabolic  model, 

Rlt^/Ri^  * - 1.  Consequently,  the  typical  value  of  the  normalized 
vacuum  field  b^  in  the  former  case  is  close  to  unity  and  nearly 
independent  of  s,  since  the  liner  spends  comparatively  little  time  at 
large  radii  see  Fig.  3 . This  is  reflected  in  the  age  variable  5. 

This  tends  to  be  much  smaller  in  the  sinusoidal  case,  but  the  sinusoidal 
trajectory  spends  a greater  age  in  the  vicinity  of  turnaround,  where 
the  driving  field  is  near  its  maximum  value. 

When  we  use  Eq.  51  t0  evaluate  the  diffused  flux  5 as  a function 
of  time  for  these  cases,  we  obtain  the  results  shown  in  Fig.  4.  We  expect 
the  sinusoidal  trajectory  to  be  more  lossy.  This  is  seen  to  be  true, 
except  at  t « 2tQ»  Diffused  flux 


i 
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leaks  back  into  the  vacuum  region  when  the  vacuum  field  drops  below 
the  values  in  the  liner.  Since  the  flux  has  not  gotten  very  far  into 
the  liner  in  the  sinusoidal  case,  it  leaks  back  out  very  quickly.  In 
contrast,  the  "older"  parabolic  trajectories  cannot  empty  all  the  flux 
back  into  the  vacuum  region,  even  though  the  driving  field  drops  away 
more  abruptly  after  turnaround.  Note  that  in  all  cases,  this  return 
diffusion  does  not  begin  until  = 0 (c.f.  Eq.  A8),  i.e.,  somewhat 


after  turnaround. 

To  display  the  dependence  of  our  results  on  the  trajectory  and 
the  vacuum  magnetic  field,  we  proceed  as  follows.  The  diffused  flux 
for  a' constant  vacuum  field  is  obtained  by  writing  bi  = 1 in  Eq.  (51): 


■vAt  B R 2 /- 

*(0)  = — Rm°/°  I de'(0-eTl/2= 


2>Jri  B R 2 


Substituting  for  Rm  from  Eq.  (27),  we  have 


$(t)  = 2 TT  < Rj2  > i/2  B (^Kt/TT)1/2, 


where 


>-{-  fit'  R^lO 


is  the  mean  square  of  the  liner  radius  as  a function  of  time.  Thus 
the  trajectory  enters  only  through  the  factor  < R3A  > 1^2.  In  the 
general  case,  when  bj.  also  depends  on  time,  either  through  Eq.  (47) 
or  otherwise,  we  can  write 
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$(t)  = 2 TT  < Ri2  >>-/2  < Bl  > (4Kt/TT)l/2, 


(59) 


where  < Bi  > is  a shape  factor,  representing  some  time  average  of  the 
driving  field.  As  a rule,  the  product  < Ri2  > x^2  <8!  > can  only  be 
obtained  by  numerical  evaluation  of  Eq.  (51)* 

To  make  contact  with  the  results  of  Knoepfel9,  we  consider  three 
cases:  (a)  the  sinusoidal  model  [Eq.  (5*01  with  e = 1,  i.e., 

Bi  = Bq  = constant  and  Ri  = Rq  = constant;  (b)  the  e = 0.001  sinusoidal 
case  described  above;  and  (c)  a static  liner  with  Ri  = Rq  but  B/Bq  = 
sin  (TTt/2tQ)  [i.e.,  the  field  which  would  result  from  letting  e -»  0 
in  Eq.  (54)].  The  corresponding  diffused  fluxes  are  plotted  against 
time  in  Fig.  5.  Curve  (a),  which  increases  without  bound,  is  just 
found  by  writing  < Ri2  > = Rq2  in  Eq.  (57).  The  result  is  the  same 
as  that  of  Knoepfel  for  the  constant  vacuum  field  case.  Curve  (b)  is 
identical  with  the  upper  trace  in  Fig.  4.  Curve  (c)  is  found  by  re- 
placing 0 with  t = t/tQ  in  (51)*  (The  difference  between  e = 0.0 

and  e = 0.001  is  undetectible  on  the  scale  of  the  plot). 

We  can  define  a flux  skin  depth  in  the  usual  fashion  as 

s<p  = *(t)/2n  RoBq.  (60) 

Substituting  $ from  Eq.  (59)  we  obtain 

s^  = y47n~  ( < Ri2  >l/2/Ro)(  < Bi  >/Bq ) 'v/ttt  (61) 

where  *y4/n  = 1.128  is  recognized  as  Knoepfel’ s9 coefficient  for  slab 

diffusion  with  a constant  driving  field.  In  a similar  fashion,  skin 
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depths  may  be  defined  for  magnetic  energy  (s^,  joule  energy  (sT), 

and  total  energy  (s  ).  In  Table  I these  skin  depths  at  t = t are 
^ o 

tabulated  and  compared  to  Knoepfel's  values  (Ref.  9,  p.77)  for  a 
sinusoidal  driving  field  on  a slab  liner. 


I 


Figure  6 shows  plots  of  the  thermal  and  magnetic  energy  of 
diffusion  associated  with  the  trajectories  of  Fig.  5.  Since  the 
energies  are  quadratic  in  B,  the  enhanced  losses  obtained  with  the 

r * 

flatter  sinusoidal  trajectories  are  more  pronounced  here. 

For  overall  energetic  consideration,  two  measures  of  the  losses 
associated  with  flux  diffusion  are  particularly  useful  (Fig.  7)*  These 
are  the  total  (W  + W^,)  energy  loss  at  turnaround,  divided  by  the  pay- 
load  energy  Wq  = irR^Bo2/2n0  and  the  thermal  energy  at  return  divided 
by  Wq.  The  first  tells  how  much  the  finite  conductivity  affects  the 
maximum  compression,  while  the  second  measures  the  tendency  of  the 
liner  not  to  rebound  exactly,  as  well  as  the  amount  of  heat  it  absorbs 
per  cycle. 

Finally,  for  some  applications  a detailed  knowledge  of  the  diffused 
flux  profile  is  needed.  Figure  8 shows  how  such  profiles  evolve  for  a 
particular  case.  Note  that  the  picture  of  a thin  layer  near  the  boundary 
is  weakest  near  turnaround.  The  variable  T],  of  course,  measures  area, 
not  displacement.  Even  though  the  trace  for  t = 2.0  is  broader  then 
the  earlier  ones,  the  thickness  of  the  layer 


r 

i 

I 


A ~ | < T|  > Rm"l/2  Rq2/Ri  , 

where  < T]  > is  an  average  or  typical  value  of  11  for  the  profile 
evidently  drops  rapidly  after  turnaround. 
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V.  Conclusion 

The  principal  result  of  this  paper  has  been  the  establishment  of 
a simple  procedure  to  handle  diffusion  problems  in  imploding  cylindrical 
systems  at  large  Reynolds  number.  The  particular  application  treated 
here,  flux  diffusion  in  imploding  cylindrical  axial  flux  compression 
devices,  has  emphasized  the  boundary  layer  character  of  this  problem. 

If  the  equations  are  cast  in  Lagrangian  variables,  the  large 
Reynolds  number  limit  leads  to  a "slab"-like  system  of  equations  if 
the  time  variable  is  suitably  weighted  over  the  implosion  trajectory. 

To  lowest  order  in  Rm  1^2}  the  magnetic  flux  payload  is  compressed 
adiabatically  (supplying  the  boundary  condition  for  the  magnetic 
diffusion  problem).  The  diffused  flux  and  magnetic  and  thermal 
energies  are  all  confined  near  the  compression  surface  within  a 
boundary  layer  whose  depth  is  of  order  Rm-1^2. 

The  main  aspects  of  this  approach  are  also  valid  in  other 
problems  involving  the  diffusion  of  a scalar  quantity  into  a moving 
medium  at  large  associated  Reynolds  numbers,  that  is,  wherever  the 
hydrodynamic  time  scale  dictated  by  the  basic  system  motion  is  short 
compared  to  the  relevant  diffusive  time  scale.  The  same  procedure 
can  treat  slab,  cylindrical,  and  spherical  diffusion  problems  in  this 
limit . 
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Appendix:  The  Convection-Diffusion  Problem  in  Slab  Geometry 

As  might  be  expected,  the  planar  analog  of  the  problem  of  flux 
diffusion  into  a cylindrical  liner  is  relatively  tractable  and  affords 
a convenient  means  of  understanding  the  principal  physical  features 
of  the  latter.  Now  the  vacuum  region  is  located  between  a perfectly 
conducting  wall  at  x = 0 and  a moving  liner-vacuum  interface  at  x = xj 
(see  Fig.  2).  Equations  (1 6 -21  ) become 


§ = (x  - x1)/xq; 


&b 

6t 


b(5,0)  = 0,  ? > 0 ; 


lim  b ( ? , t ) = 0 ; 

^ CO 


blXl/xo  + f d5  b(?,T)  - i, 

O 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 


with  xq  = xi(0)  and 

Rm  = x 2/k  t . 
o o o 


(A6) 


If  bi(T)  were  known,  the  solution  of  Eqs.  (A2-A4)  could  be 
written  down  immediately  using  standard  results  from  diffusion  theory^0 
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with  t defined  by  Eq.  (15)* 

The  Abel  equation  (A10 ) can  be  solved  by  Laplace  transforms, 
provided  iu  > 0.  Define 

00 

a(s)  =J  dcre-SCT  bi(a).  (13) 

o 


Then 


i(l  + s"l/2)  + 


da 

O)  — : = S 

ds 


(Al4) 


whence 


a(s  ' 


"t/ 


dX 


exp 


s(xg-i)+a/s  ( X— l ) 

(JO 


(A15) 


Inversion  of  (A15)  leads  to 

Vj.+uoa 

bi(a)  = 


dX  X-l 


J-n  ao2 


ttlf:  exp 


■ (X-  1 )a/(l)g 


a + 


1-X‘ 


(A  16) 


Next  we  note  (by  direct  substitution  if  necessary)  that  the  solution 
(Al6)  works  equally  well  for  uu  < 0,  provided  1 + uxr  > 0.  Finally  we 
reintroduce  t and  Rm  explicitly,  rewriting  the  integral  in  the  process 
to  obtain 


bi(T)  = 


— — f dz  e 

^^T^tr  J 


ID 


\/l+(l-Hi)T  )(z2-l)[  1 - 2/  1+  ( 1-Hut  )(z2-l )] 


(A17) 


a form  which  is  only  valid  for  u)  < 0 . 
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Though  unwieldy,  this  result  provides  some  useful  information  • 


First,  the  limit  m -»  0 yields  the  static  liner  result 


t>i(T)  = eTerfc(v/T)  (A1 

where  erfc  denotes  the  complimentary  error  function.  For  t -»  <»,  the 
asymptotic  form  of  the  latter  leads  to 


bi(T) 


(Aly 


Secondly,  when  1 + cut  = 0 (totally  imploded  liner),  Eq.  (A17)  reduces 
to 


b1(-(u“1)  = 


T 

tt 


dz  e z(l+z) 


(A2(J 


Thus  the  vacuum  field  strength  remains  finite  no  matter  how  rapid  the 
compression;  indeed,  if  Rm  is  less  than  a critical  value  Rm  given  by 


Rm  = | [l  + i (1  + Jl  + 2tt)J  «*»  1.09 


(A2: 


bi  decreases  monotonically.  Lastly,  for  ou  » 1 (the  large  magnetic 
Reynolds  number  limit), 


bl(T)  1+5F  xo/xi  ‘ 


Corrections  to  the  lowest  order  result  are  evidently  ~ Rm"1/2  and  may 
be  calculated  systematically. 

In  the  limit  Rm  » 1,  however,  Eq.  (A5)  reduces  to  b1(T)  = x /x1? 
independent  of  any  ansatz  of  the  form  (A9).  Given  a knowledge  of  the 
trajectory  x^t),  all  of  the  quantities  of  physical  interest  can  be 
found  from  (A7).  Thus  the  total  flux  diffused  into  the  liner  is 


MO 


/B  X T bi(T 


The  rate  of  change  of  $ with  time  is  given  by 


00  o 

-/ d,)  / 


rn  - o 


Hence  from  (AJ), 


i.  5b 

T|-0+  ~ 


bi(T') 


= - J-  * /V  — 

a/tt  J VT-t 


bi(0)  . -T  Sbi/bT' 

L.  ( dT  

/ttt  Jri  J (t-t')1/2 


The  magnetic  energy  associated  with  the  diffused  flux  is 
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WM  = 


x B 


2» 


b/  d,lb2(V) 

I o 


(A 26) 


B ^x 
o o 


MOMt"  ) 


TT  ^/2  /*  , , , bl(T  )bl( 

-5-  ) / dT  / dT  - =====  = 

J J0  J* 


after  an  interchange  in  orders  of  integration.  The  Joule  heat  is 


(A27) 


The  interchange  of  orders  of  integration  is  not  valid  here.  Instead  we 
integrate  by  parts,  using  (A25): 


wT  = 


x B " 
o o 


^ )^_/V  [bl(0,o 


6bx" 


(A28) 


x B „ 

O O , TT_ 

U Rm 


dT*bl(0 

o 


bbi/fcT' 

"(T*-Vr  b2 


Although  (A26)  and  (A27)  are  manifestly  positive,  (A28)  is  not. 
Nonetheless,  numerical  evaluations  of  (A2 6)  and  (A28)  reveals  that 
WT  2 W , a relation  which  can  be  shown  to  hold  in  general. 
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TABLE 

I SKIN  DEPTHS  FOR 

SINUSOIDAL 

TRAJECTORY:  Comparison  with 

Static  Liner  and 

Sinusoidal 

Field  at  t = t 

0 

Flux 

Joule 

Magnetic 

Total 

e 

s /^t 

s AAT 

T 

s A/RT 
M 

s A/^t 
E 

.001 

1.065 

.684 

.582 

1.266 

.01 

1.067 

.686 

.583 

1.270 

.15 

1.081 

• 753 

.598 

1.351 

•30 

1.092 

.856 

.611 

1.468 

•45 

1.101 

.984 

.625 

1.606 

Slab9 

.88 

• 718 

• 500 

1.223 
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INITIAL  POSITION 
OF  LINER-VACUUM 
INTERFACE 


LINER-  X 
VACUUM  \ 
INTERFACE  \ 
AT  TURNAROUND 


‘R°y 


INSTANTANEOUS 
LINER-VACUUM 
INTERFACE  AT 
TIME  T 


Fig.  1 — Schematic  of  cylindrical  liner  compressing  an  axial  flux. 
The  radius  of  the  liner-vacuum  interface  is  Rj  = Rce— 1/2  at  t = 
0 and  Rj  = R0  at  turnaround  (t  = tQ).  The  shaded  region  rep- 
resents the  layer  occupied  by  diffused  flux. 
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INFINITE  §§ 
CONDUCTIVITY 


FINITE 

CONDUCTIVITY 


VACUUM 


Fig.  2 — Slab  equivalent  of  the  model  shown  in  Fig.  1.  A region  of  in- 
finite conductivity  lies  to  the  left  of  x = 0.  Inside  the  “liner”  (the 
region  of  finite  a,  x > x3 ) the  field  strength  falls  from  its  vacuum  value 
to  zero. 
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F.  Perkins' 

D.  Grove 
H.  Okuda 
P.  Rutherford 
J.  M.  Greene 

D.  Jassby 
W.  Tang 
H.  Furth 
R.  Kulsrud 
R.  Mills 
C.  Oberman 

J.  File 
E.  Frieman 
J.  Johnson 
D.  Meade 
T.  K.  Chu 
T.  Stix 

•} 

! 
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D 


Universities 


1.  College  of  William  and  Mary 
Williamsburg,  Virginia  23185 

ATTN:  S.  P.  Gary 

Fred  Crownfield 

2.  Cornell  University 
Ithaca,  New  York  14850 

ATTN:  Prof.  R.  N.  Sudan  (App.  Phys.) 

Prof.  E.  Ott  (E.  E. ) 

R.  V.  Lovelace 
C.  Wharton 

3.  Institute  for  Advanced  Study 
Olden  Lane 

Princeton,  New  Jersey  O854O 
ATTN:  M.  Rosenbluth 

4.  Massachusetts  Institute  of  Technology 
77  Massachusetts  Avenue 
Cambridge,  Mass.  02139 

ATTN:  B.  Coppi  (Phys.)  L.  Lind sky 

J.  McCune  (A.  E. ) D.  Rose 
A.  Bers  (E.  E.)  T.  H.  Dupree  (Nuc.  Eng.) 

5.  New  York  University 

Courant  Institute  of  Math.  Sci. 

251  Mercer  Street 

New  York,  New  York  10012 

ATTN:  H.  Grad 

W.  Grossman 
H.  Weitzner 
J.  Tataronis 

6.  Stevens  Institute  of  Technology 
Hoboken,  New  Jersey  07030 


7- 


ATTN:  George  Schmidt 

B.  Rosen 

University  of  California 
Berkeley!  California  9*^720 
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8.  U.C.L.A. 

Los  Angeles,  California  9002k 

ATTN:  J.  Dawson  (Phys. ) 

Y.  C.  Lee 

9.  University  of  Maryland 
College  Park,  Maryland  20742 


B.  Fried  (Phys.) 
F.  Chan  (E.  E. ) 


ATTN:  P.  C.  Liewer  C.  S.  Liu 

Prof.  R.  Davidson  H.  Griem 


10.  University  of  Rochester 

Rochester,  New  York  14627 

ATTN : Moshe  Lubin 

L.  Goldman 
P.  J.  Gat to 


A.  Simon 
R.  McCrory 


11.  University  of  Texas 

Austin,  Texas  ’J&J12 


ATTN: 


D.  Ross 
W.  Drummond 
A.  Ware 
L.  Sloan 


F.  Hinton 
W.  Horton 
R.  Hazeltine 


12.  University  of  Wisconsin 

Madison,  Wisconsin  55706 


ATTN: 


R.  Conn 
R.  Boom 


G.  Rule inski 


15«  Yale  University 

New  Haven,  Conn.  06520 

ATTN:  B.  Bernstein 
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E.  Industrial  Research 

1.  Batelle 

Pacific  Northwest  Lab. 
p.  o.  Box  999 

Richland,  Washington  99352 
ATTN:  L.  Schmidt 

2.  Bell  Telephone  Laboratories 
Whlppany,  New  Jersey  07981 

ATTN:  S.  Buchsbaum  A.  Hasegawa 

3.  Forsen  Incorporated 
777-106  Avenue,  2*.  W. 

Bellevue,  Washington  98004 

ATTN:  H.  Forsen 

4.  General  Atomic  Company 

P.  0.  Box  81608 

San  Diego,  California  92138 

ATTN:  T.  Ohkawa  C.  Baker 

D.  Dobrott  J.  Purcell 

G.  Guest 

5.  Maxwell  Laboratories,  Incorporated 

9244  Balboa  Avenue 

San  Diego,  California  92123 

ATTN:  A.  Kolb  A.  Trivelpiece 

V.  Fargo  N.  Rostoker 

6.  Physics  International  Company 

San  Leandro,  California  94557 

ATTN:  V.  Bailey 

7.  Sandia  Laboratory 

Alburquerque,  New  Mexico  8fll5 

ATTN;  Thomas  Wright  J.  Freeman 
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8. 


Science  Applications,  Inc. 

Laboratory  for  Applied  Plasma  Physics 

Research  Staff 

La  Jolla,  California  92037 


ATTN:  N.  Krall 

R.  Sharmy 
J.  McBride 


C.  Wagner 
H.  Klein 
N.  Byrne 


9.  Westinghouse  Electric  Company 
Breeder  Reactor  Division 
P.  0.  Box  355 

Pittsburgh,  Pennsylvania  I523O 

ATTN:  Z.  Shapiro  D.  Klein 

R.  Rose  I.  Varljen 


10.  Fusion  Energy  Corporation 
P.  0.  Box  2005 
Princeton,  N.  J.  085^0 

ATTN:  B.  C.  Maglich  R.  A.  Millar 
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International 


Australian  National  University 
Canberra  A.  C.  T.  Australia 

ATTN:  J.  D.  Strachan  R.  L.  Dewar 

Center  d'Etudes  Nucleaires 
Boite  Postale  No.  6 
92260  Fontenay-Aux-Rose 
France 


ATTN:  J.  Tachon 

C.  Mercier 


N.  Troche ri a 
R.  B.  Paris 


Centro  di  Studio  sui  Gas  Ionizzati 
Unlversita  di  Padova 
6/a  Via  Gradenigo 
35100  Padova,  Italia 

ATIN:  G.  Malesani  G.  Nalesso 

Culham  Laboratory 
Abingdon,  Berkshire,  England 


ATTN: 


Chris  Lashmore  Davies  R. 


I.  Cook  M. 
K.  Roberts  J. 

J.  Christiansen  T. 
J.  B.  Taylor  J. 
B.  Peas  A. 
R.  Rebut 


Bickerton 
Hughes 
A.  Wesson 
Stringer 
W.  M.  Paul 
Gibson 


Institute  of  Plasma  Physics 
Nagoya  University 
Furo-cho,  Chikusa-ku 
Nagoya  464,  Japan 

ATTN:  S.  Hiroe  T.  Dawamura 

H.  Ike garni  T.  Kawabe 

Japan  Atomic  Energy  Research  Institute 
Tokal-mura,  N aka- gun 
Ibaraki  319-11 
Japan 


ATTN:  N.  Fujisawa 


8. 


9. 


10. 


n. 


12. 


13. 


Kurchatov  Institute 
Moscow,  U.  S.  S.  R. 


ATTN:  B.  B.  Kadomtsev 

L.  I.  Rudakov 


E.  P.  Velikov 


Laboratoria  di  Fisica 
Via  Celoria  16 
20133  Milano,  Italia 


ATTN:  P.  Caldirola 


E.  Sindoni 


Max  Planck  Institute  fur  Plasma  Physics 
8046  Garching  bei  Munchen 
West  Germany 


ATTN:  D.  Biskamp 

Horst  Pacher 
W.  Schneider 


K.  Von  Hagenow 
Chodura 
D.  Pfirsch 


Siberian  Nuclear  Physics  Institute 
Novosibirsk,  U.  S.  S.  R. 


ATTN:  D.  D.  Ryutov 


University  of  Natal 
King  George  V Avenue 
Durban,  South  Africa 


ATTN:  M.  A.  Hellberg 


CNEN,  Laboratory  Gas  Ionizzati 
C.  f.  65^  00044  Frascatt: 

Rome,  Italy 


ATTN:  R.  Toschi 


FOM-Institut  voor  Piasmafysik 
Rijnhuizen,  P.  0.  Box  7 
Jutphaas,  Netherlands 


ATTN:  C.  M.  Braams 


F.  Engeiman 


- 


9 


